Field Dependence of the Electron Spin Relaxation in Quantum Dots 
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Interaction of the electron spin with local elastic twists due to transverse phonons has been studied. 
Universal dependence of the spin relaxation rate on the strength and direction of the magnetic field 
has been obtained in terms of the electron gyromagnetic tensor and macroscopic elastic constants 
of the solid. The theory contains no unknown parameters and it can be easily tested in experiment. 
At high magnetic field it provides parameter-free lower bound on the electron spin relaxation in 
quantum dots. 
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Relaxation of the electron spin in solids is a funda- 
mental problem that has important applications. Among 
them are electron spin resonance and quantum comput- 
ing. In a semiconductor quantum dot the relaxation time 
for the electron spin is determined by its interaction with 
phonons, nuclear spins, impurities, etc. While impurities 
and nuclear spins can, in principle, be eliminated, the 
interaction with phonons cannot. Thus, spin-phonon in- 
teractions provide the most fundamental upper bound on 
the lifetime of electron spin states. The existing methods 
of computing electron spin-phonon rates in semiconduc- 
tors rely upon phenomenological models of spin-orbit in- 
teraction, see, e.g., Refs. EUlllllllBllllllE These 
models contain unknown constants that must be ob- 
tained from experiment. Meantime, as has been noticed 
more than 50 years ago by Elliot (see also Ref. H^ . 
the spin-orbit coupling in semiconductors determines the 
difference of the electron g-factor from the free electron 
value of 50 — 2.0023. The question then arises whether 
the effect of the spin-orbit coupling on spin-phonon re- 
laxation can be expressed via the difference between the 
electron gyromagnetic tensor gap (a, P = x,y, z) and the 
vacuum tensor goSap. Since gap can be measured inde- 
pendently, this would enable one to compare the com- 
puted relaxation rates with experiment without any fit- 
ting parameters. In this Letter we show that this, in- 
deed, can be done under certain reasonable simplifying 
assumptions. 

Zeeman interaction of the electron with an external 
magnetic field, B, is given by the Hamiltonian 



has been studied in the past [3|, IJ, Ig by writing all terms 
of the expansion of gap on the strain tensor, Uap^ per- 
mitted by symmetry. This gives spin-phonon interaction 
of the form Aap-^sUapo-yBs with unknown coefficients 
Aap-yS- To avoid this uncertainty we limit our considera- 
tion to local rotations generated by transverse phonons. 
The argument for doing this is three-fold. Firstly, the 
rate of the transition accompanied by the emission or 
absorption of a phonon is inversely proportional to the 
fifth power of the sound velocity dSj. Since the velocity 
of the transverse sound is always smaller than the veloc- 
ity of the longitudinal sound l^J , the transverse phonons 
must dominate the transitions. Secondly, for a dot that 
is sufficiently rigid to permit only tiny local rotations as 
a whole under an arbitrary elastic deformation, the emis- 
sion or absorbtion of a quantum of the elastic twist will 
be the only spin-phonon relaxation mode. Finally, we 
notice that interaction of the electron spin with a local 
elastic twist generated by a transverse phonon does not 
contain any unknown constants. Consequently, it gives 
parameter-free lower bound on the electron spin relax- 
ation rate. 

The angle of the local rotation of the crystal lattice in 
the presence of the deformation, u(r), is given by 



(50 = -V X U, 



(3) 



Hz — ~l^B 9a(3 SaBp , 



(1) 

(t/2 is the di- 



where /i^ is the Bohr magneton and s 
mensionless electron spin with a a being Pauli matrices. 
One can choose the axes of the coordinate system along 
the principal axes of the tensor gap- Then gap is diago- 
nal. 



9aP — gaSaP, 



(2) 



and the local angular velocity is (50. The analysis of the 
effect of the rotation on the electron spin can be done 
in the coordinate frame that is rigidly coupled to the 
crystal lattice. In that coordinate frame the effect of the 
rotation is two-fold. Firstly, it results in the opposite 
rotation of the external magnetic field felt by the spin. 
The corresponding perturbation of the magnetic field is 
given by (5B = B x (50. Secondly, the Hamiltonian in the 
rotating frame acquires a kinematic term —hs ■ (50. The 
spin-phonon interaction in the lattice frame (marked by 
prime) is then given by 



represented by three numbers, gx, gy, and gz that can be 
directly measured. Perturbation of Eq. ^ by phonons 



^'a = S^a + U^B/f^)9aP [B X (50]^ . 

(4) 
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In these formulas S(p should be understood as an oper- 
ator. Summation over repeated indices is implied. The 
canonical quantization of phonons and Eq. |(2J) yield 



2MN ^ 

kA 



2MN ^ 

kA 



OkA 

/^kA 

[ik X GkAje*^ 

/t^kA 



-kA 



(5) 



(flkA + fllkA) 



(6) 



where M is the mass of the unit cell, N is the number 
of cells in the crystal, ekA are unit polarization vectors, 
A = ti,t2,l denotes polarization, and 0Jk\ = v\k is the 
phonon frequency. The total Hamiltonian in the lattice 
frame is 



-'^s — ph ' 



Hn — Hz 



ph, 



(7) 



where Hz is Zeeman Hamiltonian of Eq. unperturbed 
by phonons and Hph is Hamiltonian of free phonons. 

Spin-phonon transitions occur between the eigenstates 
of Ho- These eigenstates are direct products of the spin 
and phonon states 



(8) 



Here \ip±) are the eigenstates of Hz with energies E± 
and \ 4>±) are the eigenstates of Hph with energies -Ephi, 
satisfying 



pii,+ 



E, 



ph, 



(9) 



For F^_ph 



of Eq. which is linear in phonon ampli- 
tudes, the states \4'±) differ by one emitted or absorbed 
phonon with a wave vector k. We will use the following 
designations 



|0+) = |n-k), \(j>_) = \nu + l) 



(10) 



We need to compute the matrix element corresponding 
to the decay of the spin |*+) |vl/_). With the help of 
Eq. igj we get: 

(Vl._|4'_pjvl;+)=K. (11) 

where components of vector K are given by 

= -f^B {9a - gp) Bptap^ \Sa\lp+) , (12) 

and the principal components of the gyromagnetic tensor, 
ga, are defined by Eq. To obtain Eq. (|12l) . we have 
used the relation 



(13) 



to eliminate dcj) from Eq. Q), and the energy conserva- 
tion, Eq. . Note that for the isotropic factor = 
and thus phonons do not couple to the spin. 



As an independent test, one can consider the problem 
in the laboratory frame. In the presence of the local 
rotation given by Eq. ^ , the gyromagnetic tensor in the 
laboratory frame becomes 



(ph) 

gap 



laa'^PP'ga'P', 



(14) 



where R^^ is a 3 x 3 rotation matrix. For small 54>, one 
has 



— Sap — ^aPjScj)^. 



(15) 



Substituting Eq. (|14|l into Eq. Q and using the orthogo- 
nality of the rotation matrix. Rap = K^^, we get for the 
Zeeman Hamiltonian in the presence of phonons 



H 



(ph) 



-Mi3.9a'/3' KL^a) (V/J^/^) ' (1^) 

In the linear order in 6<p with the help of Eq. (|15|) one 
obtains the full Hamiltonian in the laboratory frame 



H = Ho + -ffs-ph , Hq = Hz + -ffph, 



where 



and 



Hs-ph = -fi^ ■ s, 



(Mis/^) [ga - gp) Bpeapjl 



(17) 



(18) 



(19) 



One can see that the spin-phonon matrix element in the 
laboratory frame, (^-| Hs-ph |^'+) , is the same as that 
in the lattice frame, Eqs. Hll|l and H12|l . 

To obtain the relaxation rate one can use the Fermi 
golden rule. With the help of Eq. (j^Jl the transition ma- 
trix element of Eq. (|ll|l can be expressed as 

-\ Hs-ph 1"^+) 

= ^E^kAK' + llflkA + aLkAl'^k'), (20) 
kA 

where 



VkA 



:K- [k X GkA]- 



(21) 



kA 



Note that only the transverse phonons contribute to the 

relaxation process. The decay rate W ^ of the upper 

spin state into the lower state, accompanied by the emis- 
sion of a phonon, and the rate W-\ — of the inverse pro- 
cess are given by 



W+- 



(22) 



where n^^^ = (e''"o/('=sT) _ 2) 1 is the phonon occupation 
number at equilibrium. 



1/2 



\ 7 / 



(23) 
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is the distance between the two spin levels, and 

VFo = ^5]|14A|'27r5(wkA-wo)- (24) 



N 

kA 

The balance equation for normalized populations of the 
upper spin state ri+ and the lower spin state n_ (satis- 
fying ri+ + n_ = 1 ) is 

n+ = -W^+n+ + W+-n^ ^ -r n+ + W+- (25) 

where the relaxation rate is given by 

hujQ 



r = w+- + W-+ = Wo coth 



2fcsT 



(26) 



Using Eq. ||2U and replacing (1/-/V) X^k ■ ■ ■ ^ly 
vq J d^k/{2nY ... in Eq. H24|l . vq being unit-cell volume, 
one obtains 



Wo = 



1 IKI^c^i^ 

1 12 / . t3 



1 |Kpc^3 



12TThMvtuj]j 12nh pvl 



(27) 



Here vt is the velocity of the transverse sound, p is the 

1/3 

mass density, ujd = "^tVo '^^ Debye frequency for the 
transverse phonons, and |K|'^ = ^ K*K^. With the help 
of Eq. we get 

^^^lY. (ff" - 9pf {BlTo^c + BpB^T^p) , (28) 

a/3 

where 

T^p = \ ((V- kol \sp\i^+) + c. c.) . (29) 

The eigenstates of Hz are 

= -sin(6'/2)e-*'^/2|+)+cos(0/2)e*'^/2|_^ 
= cos(^?/2)e-*^/2 ^ sin(0/2)e'^/2 ^ (gg) 

where |-|-) , |— ) are the eigenstates of the operator s^. 



s.|±) = ±-|±). 



(31) 



and the spherical angles 6 and are defined through 



b = gxHxGx + gyHyBy + g^H^e^ 

= |b| (sin cos (/J -I- sin sin + cosSe^) , 



(32) 



This gives 



= — (isini/? -|- cos0cos(/9) 
("0- 1 = - (— icosyj -I- cos6'sin(^) 

{■ip-\sz \ip+) = -isin6'. 



(33) 



Direct calculation then yields 



If. gaBagiiBp 



(34) 



Thus one obtains 
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a/S—x^y.z 



{g^+gpfBlBl 



(35) 



With the help of Eq. Eq. now can be written 
in the final form 



1 {psB) 



37rft Mvl {hwof 
where n = B/_B and 



FT(n) 



^ / PbB 
Snp \ hvt 



Ft (n) , 
(36) 



i^.(n)^^(E«) 



3/2 



coth 



7 

a/3 



2kBT 



E 



1/21 



2 2 

57'^7 



(ga +ff/3)^n^ n|; 
E7 (37%) 



2 I 2 
"a +'^/3 



(37) 



Here a,/3,7 run over x,y,z. It is apparent from Eq. 
(|37|l that the relaxation mechanism studied in this Letter 
requires anisotropy of the gyromagnetic tensor. When 
the field is directed along the z-axis, Eq. (|37f) simplifies 
to 

Frie.) = |[(<?. - g.? + {g. - gyf] coth (1^) . 

(38) 

For the theory to be valid, wo of Eq. H23|l should not ex- 
ceed LJ D , otherwise there will be no acoustic phonons re- 
sponsible for the discussed spin-phonon relaxation mech- 
anism. If ga are of order unity, this is equivalent to the 
condition that the factor [p, gB / hvt) in Eq. (|36|l . that has 
dimensionality of the wave vector, is less than the Debye 
wave vector, kjj = u}£,/vt^ for transverse phonons. This 
condition is almost always satisfied in the experimentally 
accessible field range. At fuvo 3> fc^T the coth factor in 
equations 1)3 7|l and (|38|l tends to one. In this case F oc J5^ 
while Ft depends only on the direction of the field with 
respect to the principal axes of gap ■ In the opposite limit 
of fcsT ^ hivo, the relaxation rate is proportional to _B'*T 
while its dependence on the direction of the field is given 
by the factor 



E ^3o, - giif 



a/3 



2 (ga+g/3) nlnl 

E7 (57%) 



(39) 
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A nice property of the spin relaxation mechanism stud- 
ied above is its universal dependence on the strength and 
the direction of the magnetic field. Due to in Eq. H3()|) 
this mechanism can dominate electron spin relaxation at 
high fields. For, e.g., p ~ Sg/cm'^ and vt ^ 2 x lO^cm/s, 
it gives r ~ 3 X fO''s~^FT(n). The dependence of the 
rate on the direction of the field, ^^(n), is entirely de- 
termined by the difference between principal values of 
the tensor gafs- Highly anisotropic gaf3 has been theo- 
retically predicted in two-dimensional systems llfil and 
experimentally detected in GaAs quantum wells In 
the mK temperature range, the spin relaxation times of 
order 100 /iS have been observed in GaAs electron 
quantum dots in the field of order lOT. Note that equa- 
tions (|2Sl-lISni) allow a detailed comparison between the- 
ory and experiment for the proposed mechanism of relax- 
ation, which, thus, can be easily confirmed or ruled out 
for a particular quantum dot. 

In conclusion, we have studied electron spin relax- 
ation in quantum dots due to local rotations generated 
by transverse phonons. This is unavoidable relaxation 
channel that occurs when the electron gyromagnetic ten- 
sor, (7q/3, is anisotropic. It can dominate spin relaxation 
at high magnetic fields. The advantage of our theory is 
that it expresses the effect of unknown spin-orbit inter- 
actions on electron spin relaxation in terms of the tensor 
ga/s alone. The corresponding relaxation rate has univer- 
sal dependence on the strength and direction of the field 
with respect to the principal axes of gap- The important 
feature of the proposed mechanism is that it docs not 
involve any unknown constants of the quantum dot and 
is entirely determined by the three principal values of 
gaf3, which can be independently measured. This allows 
simple experimental test of the proposed theory. 
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